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A BUOYANT TORNADO- PROBE CONCEPT INCORPORATING 
AN INVERTED LIFTING DEVICE 


By Frederick C. Grant 
Langley Research Center 

SUMMARY 

Addition of an inverted lifting device is shown to make possible low-altitude entry to 
tornado cores with easier launch conditions for balloon- lifter combinations than for the 
simple balloon probe. Balloon- lifter combinations are particularly suitable for penetra- 
tion of tornadoes with average to strong circulation, but tornadoes of less than average 
circulation which are inaccessible to simple balloon probes become accessible. The 
increased launch radius which is needed for access to tornadoes over a wide range of 
circulation results in entry times of about 3 minutes. For a simple balloon probe the 
uninflated balloon must be first dropped on, or near, the track of the tornado from a safe 
distance. The increase in typical launch radius from about 0.75 kilometer to slightly over 
1.0 kilometer with a balloon- lifter combination suggests that a direct air launch may be 
feasible. 


INTRODUCTION 

Reference 1 describes in some detail a technique which makes possible the introduc- 
tion of instrumented balloons into the cores of tornadoes by means of an airplane which, 
from a safe distance, parachutes the uninflated balloon on, or near, the projected track of 
the tornado. After the tornado has closed to some 0.75 kilometer the balloon rapidly 
inflates and is driven by buoyancy forces up in the gravitational field and inward in the 
centrifugal field of the tornado. Typical times from inflation to core entry are between 2 
and 3 minutes. The self-guidance of the balloon and the safe, several-kilometer, minimum 
approach distance of the airplane are outstanding virtues of the technique. 

Dependence on buoyancy for the driving force requires that the upward, gravity- 
induced buoyancy component be accepted along with the more desired inward component. 
Thus, a finite altitude of core entry exists for every buoyancy and launch distance. That 
portion of the funnel below the' core- entry altitude is inaccessible to the probe. If the 
upward buoyancy force can be counteracted without at the same time substantially detract- 
ing from the inward buoyancy force, the entire funnel height is laid open to exploration. 
Furthermore, the analysis shows that entry times can be reduced with the lifter. 



A simple method for counteracting the upward buoyancy is to tether an inverted 
lifting device below the balloon, thereby curving the flight path downward. Many forms of 
lifter can be envisioned such as high aspect-ratio gliders, kites of low lift-drag ratio, 
parawings, and low aspect-ratio gliders such as the Rogallo paraglider (ref. 2). Yet 
another possibility is a balloon of streamline form trimmed by tail surfaces so as to be 
subject to a downward aerodynamic force. This device is the limiting case of zero tether 
length in which buoyant and aerodynamic lift are produced by a single body. Specific 
results are given for the paraglider as representative of the entire class of possibilities. 

“ The paraglider is regarded as representative only in terms of its lift and drag 
characteristics. Every form of lifter surely will have its own difficulties with stability, 
deployment, rigging, and interaction with the balloon wake. The analysis to be given pro- 
duces trajectory results valid for all forms of lifter if the ancillary difficulties are 
assumed to be overcome. 

Just as for the simple balloon, the moment of deployment of the balloon- lifter combi- 
nation is the most demanding structurally. However, the deployment of the lifter can be 
delayed by a couple of seconds, by which time the relative wind will have diminished from 
gale or hurricane levels to several meters per second. (See ref. 1.) 
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lift- drag ratio of balloon- lifter combination, /3/e 

radial and vertical components, respectively, of balloon net buoyancy force, N 

resultant net buoyancy on balloon- lifter combination, N 

tornado- strength parameter, (r/ 27 r) 2 /g, m3 

drag coefficients of balloon and lifter, respectively 

lift coefficient of lifter 

balloon diameter, m 

aerodynamic drags of balloon and lifter, respectively, N 


Dn 


,F z ,Fg 


total drag of balloon- lifter combination, N 

aerodynamic forces on balloon-lifter combination in radial, vertical, and 
tangential directions, respectively, N 
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resultant aerodynamic force on balloon- lifter combination, N 
acceleration of gravity, 9.80665 m/s2 
trajectory integral over radius, m~2 

aerodynamic lift of lifter, N 
lift-drag ratio of lifter 
mass of lifter, kg 

total inertial mass of balloon- lifter combination, kg 
virtual mass of balloon, kg 
atmospheric pressure, N/m2 
ratio of balloon drag to lifter drag 

ratio of balloon cross-sectional area to lifter planform area 

radial distance from tornado center, m 

radius at peak altitude of trajectory, m 

launch radius, m 

frontal area of balloon, m2 

planform area of lifter, m2 

tether tension, N 

time to center 

wind speed, m/s 
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r 
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. Ap 
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i Ap T 


balloon volume, m3 

velocity in radial, vertical, and tangential directions, m/s 

relative wind, m/s 

gravitational force on lifter, 

centrifugal force on lifter, N 

altitude, m 

altitude of probe at center, m 
altitude of probe at trajectory peak, m 
density scale height, 8436.1 m 

i2 

tornado circulation, m^/s 
circulation of mean tornado, 1.75 x 10^ m^/s 
flight-path angle with horizontal, radians 
balloon mass density of buoyancy, kg/m3 
total mass density of buoyancy, kg/m3 


e = C D,K S K + C D,B S B» m2 
8 azimuthal coordinate, radians 

X = arc tan a, radians 
u = arc tan — , radians 

' rrr » 7 


l,S 
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radial and vertical components, respectively, of square of relative 
wind, m2/s2 



v = arc tan (L/D)j£, radians 


P 

P B 

% 

PT 

p o 

* 


atmospheric density, kg/m3 

mean mass density of balloon, kg/m3 

equivalent density of lifter mass over balloon volume, kg/m3 
density of total probe mass spread over balloon volume, kg/m 3 
atmospheric mass density at sea level, 1.225 kg/m3 
tether angle with vertical, radians 


A dot over a symbol indicates the derivative of the quantity with respect to time. 


ANALYSIS 


Definition of Model 

The dynamical behavior of balloon- lifter combinations is investigated for the same 
simple, vertical, line-vortex tornado in an isothermal atmosphere that was used in refer- 
ence 1. The lifter adds five degrees of freedom to the system, two for the tether and three 
for the lifter, but these extra coordinates play no role in the present analysis because the 
lifter is assumed to ride at constant angle of attack with no sideslip, yaw, or roll. The 
plane of symmetry of the lifter thus always coincides with a plane through the balloon and 
the line vortex, and the lift and drag coefficients remain constant. The tether is assumed 
to lie in the plane of symmetry of the lifter and the flight-path angle is assumed the same 
for both balloon and lifter so that both degrees of freedom of the tether are eliminated. A 
view of a balloon in combination with the paraglider in flight is shown in figure 1. If the 
harness is so arranged that only small excursions in the intersection of the projected tie 
on. the chord can occur, then only small changes in lift and drag coefficients can occur. 

Lateral and longitudinal oscillations are neglected in the analysis because they have 
but a slight effect on the mean behavior. In practical cases, a lifter of low (L/D)j^ of, 
say, 4 and a sufficiently long tether between balloon and lifter should effectively inhibit 
the growth of small oscillations. (See ref. 3.) 

Finally, in addition to numerical results, a closed form, quasi-equilibrium trajec- 
tory solution such as that of reference I is found. Very quickly after launch (several sec- 
onds) certain of the inertia terms of the motion equations become small compared with the 
other terms and may be ignored with little error. The approximation is similar to that 
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used in calculation of parachute trajectories (ref. 4). For parachutes, equilibrium, or ter- 
minal, speed is often used at every air density (altitude) instead of the very slightly higher 
actual airspeed. Such a quasi-equilibrium approximation has the great advantage of reduc- 
ing the second-order motion equations to first-order equations. 

As in reference 1, the tornado winds are approximated outside the core by a verti- 
cal line vortex which produces a wind speed v equal to r/27rr for a tornado circulation r 
at a radius r. , The atmospheric density is approximated by 

p ' p o exp (- %) 

for ground-level density and density scale height Zp. The corresponding vertical pres- 
sure gradient is taken as given by the barometric equation 


where g is the acceleration of gravity. Because the Coriolis force is small compared 
with the centrifugal force, cyclostrophic balance is assumed in the radial direction; hence, 
the radial pressure gradient is 

9p _ pv2 

8r r 

The altitudes encountered during probe ascent are small enough that the exponential 
density approximation is adequate. The singularity at r = 0 is not a good approximation to 
the actual atmosphere, but use of the simple, singular vortex model does not produce sig- 
nificant changes in the trajectories or times to center. 

Equations of Motion 

The complete inhibition of lateral and longitudinal oscillations is tantamount to elim- 
ination of the several seconds of transient motion just after launch during which the probe 
rapidly catches up to, and then drifts with, the local tornado wind. In other words, for the 
tangential motion 

re = v Q = v _ (1) 

where v fl is the probe's tangential velocity. Equation (1) is the first-order approximate 

U 

equation of motion for the azimuth angle 0. The first-order approximate equations for 
the radial velocity, r = v r , and the vertical velocity, z = v z , are found by suppressing the 
v r ,v z terms in the radial and vertical equations of motion. 

The more exact motion equations of reference 1 included the virtual, or additional 
apparent, mass of the spherical balloon. When the lifter is present, corresponding equa- 
tions may be written as 
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-(m T + m')v z = -pgV B + m T g + F z (2a) 



where m-j. is the total mass of the balloon-lifter combination, m' is the virtual mass of the 
spherical balloon (half the mass of displaced fluid), and V B is balloon volume. Only the 
apparent mass of the balloon appears in equations (2). For simplicity, the components of 
apparent mass of the lifter are omitted. The regime of validity of equations (2) is the low- 
acceleration period after the lifter is deployed, when the apparent mass of the balloon is 
assumed to dominate that of the lifter. 

As mentioned previously, Vq approximates v very closely and very quickly after 
launch. Because both sides of equation (2c) vanish for vq = v, approximate equations are 

(m T + m')v z = Ap T gV B - F z (3a) 

I 2 

(m T + m’)v r = -A p t ~~ V b - F r ^ (3b) 

In equations (3) the total inertial mass of the balloon-lifter combination has been converted 
into an artificial density p^ defined by 

P T V B 

with 

A P T = P - P T 

indicating the net density of buoyancy. Because the probe speed is up to the speed of the 

2 

local wind, the additional inward force -m’ discussed in references 1, 5, and 6 cancels 

v 0 2 

against a corresponding inertial term -m’ The final step in obtaining the quasi- 

equilibrium motions is to set to zero the left-hand sides of equations (3) which are small 
compared with terms of the right-hand sides. The forces on the spherical balloon are 
indicated in sketch A. 
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Sketch A 

For motion up and to the left 
B r = (p - P b )t- V B 

and 

= (p - P B ) gV B 


are net radial and vertical buoyancies including the radial and vertical forces on the bal- 
loon in the centrifugal and gravitational fields. In sketch A, T is the tether tension, Dg is 
the balloon drag, y is the flight-path angle, and ^ is the tether angle. 



Sketch B 


In sketch B, Lx and Dx are the lift and drag on the lifter, respectively; W z x and 
w r,K are the forces on the lifter in the gravity and centrifugal fields, respectively. 

The quasi- equilibrium equations of motion may now be written as the static equilib- 
rium of the forces indicated in sketches A and B. For the vertical direction, 
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B z - T cos if/ - D b sin y = 0 

(Balloon) 

(4 a) 

T cos if/ - W Z) j£ - Dk sin y - Lk cos y = 0 
and for the radial direction, 

(Lifter) 

(4b) 

Dg cos y - T sin if/ - B r = 0 

(Balloon) 

(5a) 

T sin if/ + Dj£ cos y - Lj£ sin y + W r ^ = 0 

(Lifter) 

(5b) 


By addition in each pair of equations (4) and (5), the tether tension and tether angle can 
be eliminated. Further simplification can be made by introducing the quantities £ and £ 
where 


I = -vj el cos y ( 6a ) 


C " v rel sin r 
v rel = (* 2 + ? 2 ) 1/4 


(6b) 

(6c) 


and v re i is the relative wind. The relative wind appears in the aerodynamic drag and lift 
equations: 


D B = \ P v rel S B C D,B , 

! 

D K = 2 P v ?el S K C D,K ' 

I 

L K = ^ P v ? e l S K C L,K • 

* r ' j 

The balloon is assumed to be spherical so that the frontal area is 


(7a) 

(7b) 

(7c) 


S B =fd2 


where d is the balloon diameter. The volume of the balloon is 


V B = | dS B 


in terms of the frontal area. No special assumptions are made for the configuration of the 
" lifter; however, it is convenient to introduce into the equations of motion a fictitious den- 
sity associated with the lifter defined by 


- 3 W z,K 


PK~ 2 


gdS B 


This equation represents the density which the mass of the lifter would have if it were 
imagined to be spread uniformly over a volume equal to that of the balloon. 
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Introduction of the variables £ and £ and the aerodynamic coefficients C B g, C Bj k, 
and Cg k f rom equations (7) into the equations found by addition of equation (4a) to equa- 
tion (4b) and equation (5a) to equation (5b) yields the following pair of equations: 

-Cl,kSk£ + (c D) k s k + c D)B s B )s = | Ap T f S B (8a) 

( C D,K S K + c D,B S b)^ + C L,K S K£ = a Pt T~ £ S B ( 8b ) 

in which the relative wind v re j and the flight-path angle y do not appear. The quantity Aprp 
is defined by 


A P'f ~ P " ^P B + Pk) “ P " Pt 


where p B is the mean density of the balloon. The altitude at which Ap T = 0 is called the 
float altitude of the balloon. Solving equations (8) for £ and £ yields 




Aprp 


(9a) 


2 

, _ 4 s d ~~F C L,K 8 K + g(CD, K S K + C D,B S B) 
3 P ( C D,K S K + c D,B s b) 2 + C L,K S K 2 


(9b) 


At any radius r and density p, equations (9) yield the quasi-equilibrium velocities and 
flight-path angle through 


v rel 


= (£ + C 2 ) 


1/4 


r = v T 


z = v z = 


v rel 

C 


v rel 
Y = tan _1 (-^/ £) 


(10a) 

(10b) 

(10c) 

(10d) 


The weight of the lifter appears in Ap,p of equations (9) as p^, an addition to the balloon 
density pg. The balloon- lifter combination moves along its trajectory as though all the 
weight were in the balloon. Only the tether angle and tether tension are affected by the 
apportionment of weight between balloon and lifter. 

-The tether angle can be found by -equating to- zero the moments of-the forces on the 
balloon- lifter system. Choice of moment center at the lifter yields 
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tan i// = 


(11a) 


2 + C D,B S B^ 

2 -p^ S " C D,B S b£ 
for 

Ap B = p - p B 

The tether tension may be found from equation (4a) as 


(Hb) 


T = | pgdS B (^E - | ^ ^sec * (12) 

The factor to the left of the parentheses on the right-hand side of equation (12) is the 
weight of air displaced by the balloon; hence, the rest of the right-hand side is the tether 
tension in units of that weight. Equations (1), (10b), and (10c) are those equations which 
are numerically integrated to find the complete trajectories and entry times. 


Closed-Form Trajectories 

In the limit of no lifter = Wj- = oj, equations (9) reduce to 


t = 4 Ap B v 2 d 

5 O r* 


3 P r Cd,B 
d 


3 P ^D,B 
The ratio of equation (13b) to (13a) is 

£ = z = dz = _ gr = -gfM' 2 


(13a) 

(13b) 


or 


4 r dr v 2 


_ r 3 = 1/lr bz 


(14a) 


g\27r/ dr 


(14b) 


Equation (14b) is the balloon trajectory equation of reference 1. Similarly, the ratio of 
equation (9b) to (9a) yields 


. bz 

dr 


- -jr 0 + ge 

g/3 +j L t 


(15a) 


where 


0 = C L,K S K 

e = Cd,kSk + C D,bSb 
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Independence of the trajectories from Ap,p shows, as in the nonlifting condition, that only 
the time of traverse depends on Ap,p to which the driving force (buoyancy) is proportional. 
Introduction of the lift-drag ratio of the balloon-lifter combination 


a - S. - Ci - kSk 

6 C D,K S K + C D,B S B 

and the tornado- strength parameter 


(15b) 



where v = T/2'nr has been used, yields 



I 


dr 

b + ar3 



r^dr 
b + ar3 


(15c) 


(15d) 


where the parameter a is a measure of the effectiveness of the balloon- lifter combination 
in curving downward the trajectory and b is a measure of the relative strengths of the 
radial and vertical buoyancies. Specifically, b is the ratio, at unit distance from the vor- 
tex, of radial to vertical buoyancies. The integrals on the right-hand side of equation (15d) 
are elementary and quadrature yields 


for 


where 



I(r) 




(r + k) 2 
- kr + k 2 



k = 


/b\V3 _ 

i / r \2 C D,K S K + C D,B S B 

\a/ 

g\2ir) Cl,K s K 


1/3 


(16a) 


(16b) 


(16c) 


Altitude Peak 


Equation (15a) yields a simple result when the left-hand side is set to zero. A peak 
is found in the trajectory (as viewed in the moving plane of the balloon and line vortex) 
when 


r p = (ab) 1 /3 = 


i(r\ 2 C L,K S K 
gVnJ C D)K S K + C D B S B 


“11/3 


(17) 
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The trajectories with lift have the form shown by the solid line in sketch C. 



Sketch C 


As the balloon without lift rises along the dotted line, its motion is first dominated 
by the nearly constant vertical buoyancy and then by the radial buoyancy which increases 
with the diminishing radius. Peak altitude is at the center, r = 0. As expected, addition 
of downward lift impedes the upward motion by curving the trajectory downward. At the 
same time the gradually increasing relative wind caused by the increasing radial buoyancy 
produces lift which increases as the square of the relative wind. The lift becomes so 
strong that eventually the flight path points downward all the way to the center. Inside of 
Tp both radial buoyancy and lift are important. 

Although the balloon is descending after passing point p, it can be easily shown that 
the pressure is a monotonic decreasing function of time 

P = pgr 

Thus, the balloon remains taut throughout the trajectory and maintains constant volume. 



Force Balance on Balloon- Lifter Combination 


From equations (9) and (10a), the relative wind v re j after algebraic manipulations 
becomes 


v rel I t 


Segd 


A Pn 


A + 4\ I/4 


3 C D,K S K + C D,B S B p 


1 + a^ 


(18) 
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Equation (18) can be regarded as a statement of the balance of the total aerodynamic 
force F-p and total net buoyancy Bq\ By substitution in the numerator and denominator 
under the square root in equation (18) with the net vertical buoyancy B z and total drag of 
the balloon -lifter combination D-p, equation (18) becomes on squaring and rearranging 

(l + a2) 1 / 2 D T = B z ^l + ^) V2 (19) 

The half -power factors in equation (19) are the secants of the angles between total drag 
and total aerodynamic force and between vertical buoyancy and total buoyancy. Thus, 
equation (19) can be written 

Frp = Brp (20) 

Equation (20) and the definitions of lift and drag allow the immediate construction of the 
vector diagram of the external (no tether tension) forces on the balloon-lifter combination 
shown in figure 2. 

The motion equations (10b) and (10c) may be rewritten, by use of equations (9), (15b), 
and (15c) as 



( 21 ) 



( 22 ) 


As may be seen in equation (22), the peak altitude occurs at radius r p = ^b of the trajec- 
tory inside of which the lift and radial buoyancy dominate the motion. From equation (21) 
the corresponding r is seen to be -(v re ^p as expected. 

Squaring of equation (18) and substitution of r p = (ab)l/3 shows the dynamic pressure 
at the peak altitude to be 


(l 2 \ _ 2 s Bg d /. \ 1_ 

\ 2 PVre Vp 3 C D,B S B +C D,K S K* Pt ^P a 

which for a given lifter is independent of tornado strength b except through Ap T . However, 
Aprp varies only slightly for practical cases in which the changes in peak altitude over a 
range of tornado strength are a small fraction of the density scale height and of the float 
altitude. 
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RESULTS 


Limitations on Lift Parameter a 

The motion of the balloon-lifter combination depends on a parameter a additional to 
the tornado-strength parameter b. The form of the parameter a is 

n (V p >K 

1 + Rd 

where 


Rd 


C D,B S B 

C D,K S K 


In addition to the aerodynamic efficiency of the lifter (lift-drag ratio (L/D)^, the ratio Rd 
of ballistic parameters appears. 


In most of the trajectory calculations to follow, representative values of the aerody- 
namic coefficients are assumed. That is, Cd ( B) (L/D)k, and Cd ; k are held constant and 
RS = Sb/Sk is varied. The values of a for arbitrary (L/D)k and Rd are plotted in figure 3. 

Although the range of a in figure 3 is arbitrary, an upper limit to a must be respected 
at every launch radius. For increasingly large lifting capability the flight path will curve 
increasingly downward and for sufficiently large a-values the probe will crash short of the 
center. The limiting case is the one for which the probe crashes exactly at the center. By 
iteration in the trajectory equations (16), the limiting values of a have been calculated for 
a large number of cases and the results are plotted in figure 4. The three curves corre- 
spond to weak, average, and strong tornadoes (ref. 7). It may be fairly concluded from 
figure 4 that for launches inside of a kilometer an a-value no greater than about 1 is 
needed for strong and average tornadoes, and an a-value no greater than about 2 is needed 
for weak tornadoes. Such small values correspond to weak lifting devices. It should be 
remembered that for given aerodynamic qualities of a lifter, the value of a can be reduced 
simply by reducing the size of the lifter (increasing Rg = Sb/Sk). 

The upper values of a given in figure 4 indicate another limit on the performance of 
the balloon-lifter combination, the minimum volume in which the probe can operate from 
a fixed launch radius. This limit is defined by the altitude of the peak of the arch in the 
probe-vortex plane for the limiting a-values of figure 4. The position of the peaks for the 
three tornado strengths of figure 4 are shown in figure 5. For launches within a kilometer 
radius, probes making zero-altitude entry do not rise over about 0.3 km at any tornado 
strength. The containment power of even weak lifters is thus demonstrated. 
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• Fixed- Geometry Probes 


The results shown in figures 4 and 5 define interesting limits, but every point on 
every curve corresponds to a different value of a, hence to a different balloon- lifter com- 
bination. It is natural to consider now the performance of specific vehicles at different 
launch radii over the range of tornado strengths. 

The balloon considered is the same 2 -meter-diameter, hydrogen-filled balloon con- 
sidered in reference 1 with Cd ; B = 0-4. The lifter might be the Rogallo parawing whose 
aerodynamic characteristics may be found in reference 2. The parawing is but one exam- 
ple of many aerodynamically equivalent configurations. At high angles of attack the 
(L/D)j£ of a parawing approaches zero as it does for any wing. As the angle of attack is 
reduced, the lift coefficient rises to a value of about 0.4 at a drag coefficient of 0.1. 
Although a rigid Rogallo configuration attains higher (L/D)k values at lower angles of 
attack, a fabric configuration begins to luff (flap at the trailing edge). Values correspond- 
ing to the onset of luffing will be used. In addition to the hydrogen mass of the balloon, a 
mass of 1.6 kilograms is assumed in the lifter. Results for other masses are easily 
scaled. In practice, helium is the gas of choice; hydrogen is used in the calculations to 
define the ultimate probe performance. Use of helium results in a small decrease of 
buoyancy, 7.4 percent. 

Vehicles of different lifting capability can be studied by using different values of 
R S = Sb/Sk- An infinite value of Rg corresponds to a vanishingly small lifter. Finite 
values of Rg corresponds to different relative sizes of balloon and lifter. The largest 
lifter to be considered is one with the same area as the balloon, Rg = 1.0. The calculated 
heights at center for launch radii between 0.6 and 1.0 km in weak, average, and strong 
tornadoes are shown in figure 6. 

Nearly all of the compromises involved in choosing a size for the lifter can be dis- 
cussed in terms of figure 6. If it is desired to have the probe arrive at the core under 
the lowest cloud base, it must arrive at an altitude of no more than 0.3 km. For average 
conditions (fig. 6(b)) the base is at about 0.95 km, while the highest base is off the plot at 
about 1.5 km. For a weak tornado Rg » 2 is required for entry at average base height 
when launch is at a radius of 0.8 km. Correspondingly, Rg « 16 suffices for an average 
tornado and even a balloon alone (Rg = °°) will enter below the average base height in a 
strong tornado. Under a more stringent requirement that entry be made below the low- 
est cloud base, not even Rg = 1 suffices for a weak tornado while figure 6(c) shows that a 
balloon-lifter combination with Rg = 1 will crash in a strong tornado when launched at 
0.8 km. The conflict might be stated as: Large lifters and small launch radii for weak 
tornadoes as opposed to small lifters and large launch radii for strong tornadoes. 

A better approach might well be to abandon the weak tornadoes and to concentrate 
on low-altitude entry to average and strong tornadoes. Certainly, average and strong 
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tornadoes are more suitable subjects of study because of their more destructive charac- 
ter. Accordingly, selection of a moderate value, say Rg = 4, at a launch radius of 0.8 km 
might be considered. Such a probe would provide core entry at 0.2 to 0.6 km for average 
to strong tornadoes. The trajectories corresponding to Rg = 4 and launch at 0.8 km as 
observed in the moving probe-vortex plane are shown in figure 7. Viewed from above in 
ground-fixed coordinates the paths of the probes are always smooth spirals of constantly 
diminishing radius. Along the spirals, in the model assumed, the velocity perpendicular 
to the vortex-probe radius is always equal to that of the wind. 

The powerful effect of the lifter is evident in figure 7. Without lift in a weak tor- 
nado, the balloon alone will not ever reach the center because it comes up to its float 
altitude of 9 km first. With lift, it reaches the center in less than 6 minutes some 0.2 km 
above the maximum cloud base of 1.5 km. The time is longer and the arrival higher but 
is a great improvement nonetheless. The performance in the average tornado improves 
from arrival without the lifter at 1.3 km altitude in 3 minutes 25 seconds to arrival with 
the lifter at 0.6 km in 2 minutes 20 seconds. Less marked improvement occurs in the 
strong tornado, but the performance in a strong tornado is sufficiently impressive with- 
out a lifter. 

Launch Conditions 

At this point another influence on the choice of vehicle may be introduced; namely, 
the wind speed at launch. For the three cases presented in figure 7 which span the range 
of tornado strength, the tornado winds at a radius of 0.8 km are 12, 35, and 58 meters per 
second, respectively. The winds are thus hurricane force for the average tornado and 
about 209 km per hour (130 mph) for the strong tornado. Without experiments one can- 
not be dogmatic. However, although rapid inflation can very likely be brought off suc- 
cessfully in hurricane winds (ref. 8), it seems likely that ground launches in 209-km- 
per-hour (130-mph) winds will prove very much more difficult. Dynamic pressures are 
nearly three times as high in the strong tornado as in the average tornado at 0.8-km 
launch radius. Air launches of balloons have been successful at speeds over 161 km per 
hour (100 mph) (ref. 9). 

K sufficiently favorable launch conditions are maintained while at the same time 
a low entry height at the core is provided, the necessary cha.nge in parameters is plain: 
larger lifters and larger launch radii are needed. Consider what can be regarded as a 
shift to a very large launch radius and a very large lifter: tq =. 1.1 km and Rg = 1. The 
trajectories corresponding to these values are indicated in figure 8 in the same manner 
as are those in figure 7. 

With such a large lifter the center is reached below the mean cloud base height 
even in the weak tornado. In the strong tornado the probe crashes just short of the center. 
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Launch winds are 31, 90, and 150 km per hour (19, 56, and 94 mph) for weak, average, 
and strong tornadoes, respectively. 

Release of Lifter 

The crash indicated in figure 8 is a reminder that the core entry portion has been 
heretofore somewhat neglected as compared with the rest of the entry trajectory. At 
some point the lifter must be cut loose from the balloon so that a better response to the 
central winds of the tornado can be achieved. Although a somewhat arbitrary choice, the 
point of highest altitude seems appropriate. It marks the place where the lifter over- 
comes the vertical buoyancy and is near, but not in, the core. It happens that the peak 
altitude occurs at very nearly the same value of the dynamic pressure in the model tor- 
nadoes considered so that, for example, a flat plate drag sensor on the lifter will provide 
self-timing for cutting of the tether. The change in the trajectories resulting from cutting 
of the tether at peak altitude is indicated in figure 8. Because the probe is near the center, 
the times are not much affected (shortened slightly) but the steep descents (and the crash) 
are eliminated. 


Times to Center 

The times to center for a wide range of lifter size (the largest with Rs = 1 and the 
smallest with Rg = 16) are shown in figure 9. Some features in figure 9 should first be 
mentioned. In the weak tornado the probe with the smallest lifter is not captured in less 
than 10 minutes for radii slightly greater than 0.7 km, which is why only a single point is 
shown. The other extreme behavior, a crash short of center, first occurs in the average 
tornado for the large lifter at a launch radius of about 0.8 km, which is why the Rs = 1 
curve is cut short. In the strong tornado, the large lifter crashes for all radii inside 
about 1.1 km; hence, only a single point is shown for Rs = 1. In general, the times to 
center increase rapidly with launch radius, particularly for the smaller lifters. As with 
the desire for low-altitude entry, the desire for fast entry to the core conflicts with the 
need for acceptable launch conditions. Furthermore, if entry times are not to be far 
from, say, 3 minutes for a range from slightly below average tornado strength to the 
strongest tornado strength, then the largest lifters should be used. 

The times given in figure 9 are for a relatively light probe (a p T = 0.86 kg/m 3 ) , but 
the times in figure 9 can be easily scaled. The relative wind varies as sjAp^. Since the 
probe traverses the same trajectory if only the weight is changed (eqs. (15)), the time of 
passage varies as 
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CONCLUDING REMARKS 


Addition of an inverted, rapidly deployed, lifting device to the rapidly inflated bal- 
loon tornado probe makes possible low-altitude entry to tornado cores with easier launch 
conditions than the balloon alone. Balloon- lifter combinations are particularly suited to 
penetration of tornadoes with average to strong circulations. Although for the same launch 
conditions the time to center is shortened by the lifter, compromises made to provide 
easier launch conditions and access to tornadoes over a wide range of tornado strength 
leave the entry times in the neighborhood of 3 minutes. On the other hand, these same 
compromises increase the launch radius from about 0.75 kilometer to about a kilometer 
and suggest that an air launch of the balloon-lifter probe may be feasible. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 31, 1973. 
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Figure 3. - Variation of a of balloon -lifter tornado probe as a function of lifter (L/D)j^. 
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(a) Weak tornado, r/3. (b) Average tornado, r. 


(e) Strong tornado, 5 r/3. 


Figure 6.- Height at center as a function of launch radius. Cl k = pD B = 0-4; 
(L/D) k = 4.0; d = 2.0 m; m K = 1.6 kg; p = 0.0845 kg/m^; T = 1.75 x 10^ mVs. 














Figure 8.- Trajectories in moving plane of balloon-lifter combination and vortex 
for a large lifter at a large launch radius. Rg = 1.0; = 0.4; 

(L/D)^ = 4.0; d = 2.0 m; m*- = 1.6 kg; p D = 0.0845 kg/m3. 
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